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Blow up of smooth solutions to the barotropic compressible 
magnetohydrodynamic equations with finite mass and energy 

Olga Rozanova 

Abstract. We prove that the smooth solutions to the Cauchy problem for 
the three-dimensional compressible barotropic magnetohydrodynamic equa- 
tions with conserved total mass and finite total energy lose the initial smooth- 
ness within a finite time. Further, we show that the same result holds for the 
solution to the Cauchy problem for the multidimensional compressible Navier- 
Stokes system. Moreover, for the solution with a finite momentum of inertia 
we get the two-sided estimates of different components of total energy. 



1. The finite time blow up result 

The set of equations which describe compressible viscous magnetohydrodynam- 
ics are a combination of the compressible Navier-Stokes equations of fluid dynamics 
and Maxwells equations of electromagnetism. We consider the system of partial 
differential equations for the three-dimensional viscous compressible magnetohy- 
drodynamic flows in the Eulerian coordinates [1 for barotropic case: 

9tP + div,(H =0, (1.1) 

dt{pu) + Divx{pu (g)u) + VxP{p) = (curlj^iJ) x H + Div^T, (1.2) 

dtH — c\n:\x{u x H) = — curl2;(i^curl2:i7), divxH = 0, (1-3) 

where p, u — (1*1,1/2,^3), p, H ~ {Hi, H2, H^), denote the density, velocity, pres- 
sure and magnetic field. We denote Div and div the divergency of tensor and vector, 
respectively. Here T is the stress tensor given by the Newton law 

T = Tij = p [diUj + djUi) -\- \dvvu5ij, (1-4) 

where the constants p and A are the coefficient of viscosity and the second coefficient 
of viscosity, 1/ > is the coefficient of diffusion of the magnetic field. We assume 
that ^ > 0, A -f |/i > 0. 

The state equation has the form 

p^Ap''. (1.5) 
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Here A > 0, 7 > 1 are constants. 

Although the electric field E does not appear in the MHD system (1.1)-(1.3), 
it is indeed induced according to the relation 

E — v cnrlxH ~ u x H 

by the moving conductive flow in the magnetic field. 

The MGD system (1.1) - (1.5) is supplemented with the initial data 

{p.u,H) ^{poix),uo{x),HQ{x))en"\R''),m>3. (1.6) 

The local in time existence to the above Cauchy problem for the density sepa- 
rated from zero follows from 3 , however the results on global solvability are very 
scanty. For example, even for the one-dimensional case, the global existence of 
classical solutions to the full perfect MHD equations with large data remains un- 
solved (the smooth global solution near the constant state in one-dimensional case 
is investigated in [2]). The existence of global weak solutions was proved recently 
[4] (see also [5] for references). 



The system (1.1)- (1.3) is the differential form of balance laws for the material 
volume 0(i); it expresses conservation of mass 



m — J pdx, 
n(t) 

balance of momentum 

P = J P''^ dx, 
n(t) 

and balance of total energy 



n(t) 

Here Ek{t), E„i{t) and Ei{t) are the kinetic, magnetic and internal components of 
energy, respectively. 

If we regard fl{t) = R^, the conservation of mass, momentum and non-increasing 
of total energy takes place provided the components of solution decrease at infinity 
sufficiently quickly. 



Our main result is the following: 

Theorem 1.1. Assume 7 > |. // the momentum P ^ 0, then there exists no 
global in time classical solution to (1-1) - (1-5), (1.6) with conserved mass and 
finite total energy. 

li H = 0, we get the barotropic system of the compressible Navier-Stokes 
equations. We can consider it in any space dimensions without physical restriction 
n = 3. In this case the viscosity coefficients should satisfy the inequalities p > 
0, X + ^p > 0. The following theorem is an extension of the result of 'W to the 
larger class of solutions. 
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Theorem 1.2. Assume n > 3 and'y > -^^2- V^he momentum P ^ 0, then there 
exists no global in time classical solution to the barotropic system of the compressible 
Navier-Stokes equations with conserved mass and finite total energy. 

The next lemma is proved in [6 for both barotropic and non-barotropic cases. 

Lemma 1.3. Let n > 3, j > and u e n^{W'). If \P\ ^ 0, then there exists 
a positive constant K such that for the solutions with conserved mass and finite 
total energy the following inequality holds: 

\Dufdx> KE^""-'-'^{t). (1.7) 



In the dimension n = 3 it has the form: 

Lemma 1.4. Let 7 > |, and u e H^(R"). // |P| 7^ 0, then there exists a positive 
constant K such that for classical solutions to (1-1)- (1.5) with conserved mass and 
total finite energy the following inequality holds: 

J \Du\^dx> KE7^^{t). (1.8) 
Proof of Lemma I.4 First of all we use the Holder inequality to get 



IPI - 



pu dx 




< \ p-^dx\ \ \u\'=dx . (1.9) 



Further, using the Jensen inequality we have for 6(7 — 1) > 1 (or 7 > |) 

/ p^dx 




< 



{l-l)E.{t) 
mA 



Thus, the latter inequality and (1.9) give 

\P\ < Ki {E,{t))W^ f l^l^dx I , (1.10) 



with the positive constant Ki that depends on 7 and m. Further, we take into 
account the inequality 




<K2 / \Du\^dx, (1.11) 



where the constant K2 > depends on n, n > 3. The latter inequality holds for 
u € i?^(M") ([7], P-22) and follows from the Sobolev embedding. 

Thus, the lemma statement follows from (1.10), (1.11), with the constant K — 

|P|" 
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Proof of Theorem 1.1. Let us note that the total energy £ satisfies the mequahty 
£'{t) <-'^ J Icurl^Hp dx-cr J \Du\^ dx < 0, (1.12) 

with some positive constant a, therefore £{t) < £{Q). 
From (1.8) we have 

£'{t) < -aK{E,{t))-^i^) < -aK{£{0)y^»^ = con.3t < 0, (1.13) 

this contradicts to the non- negativity of £{t). Thus, Theorem 1.1 is proved. 

Proof of Theorem 1.2 is similar. Estimate (1.7) holds in any dimensions greater 
or equal than 3, therefore we can apply Lemma 1.3 to the case _ff = to obtain 

£'{t) < ~(jK {E,{t)y^^T^ < -aK {£{0)y^^TT^ = const < 0. 

2. Estimates of energy for solutions with finite moment of inertia 

In this section we are going to obtain more specified estimate for the decay rate 
of total energy for more restricted classes of solutions. 
Let us introduce the functionals 

Git) = \ J p{t,x)\x\'^ dx, F{t) = J {u,x)pdx, 

R3 R3 

where the first one is the momentum of inertia, the inner product of vectors is 
denoted as (., .) . 

Taking into account the general Stokes formula one can readily calculate that 
for classical compactly supported solution to (1.1)-(1.5) the following identities 
hold: 

G'{t)=F{t), (2.1) 

F'it) = 2Ekit) + Era{t) + 3(7 ~ l)E,{t). (2.2) 

In fact, formulas are obtained in [5] for gas dynamics equations with finite moment 
of inertia (not especially compactly supported). However, evidently that (2.1), (2.2) 
are satisfied for solutions decaying at infinity sufficiently quickly together with their 
derivatives. One can impose a special rate of decay on the solution and then prove 

(2.1) , (2.2) as a lemma (as in [6]). The alternative way is to take (2.1), (2.2) as a 
definition of the rate of decay. 

Definition 2.1. We will call a solution {p,u,H) to the Cauchy problem (1.1) 
- (1.5), (1.6) highly decreasing at infinity, if p{t,x) > and the identities (2.1), 

(2.2) hold. 

Evidently, the integrals of mass, momentum and total energy are finite for 
highly decreasing at infinity solutions at any fixed moment of time. 

We firstly get two-sided estimates of G{t). 

Lemma 2.2. For classical highly decreasing at infinity solution to (1.1) ~ (1-5) 
the following estimates hold: 
if J < I , then 

— t"^ + F{0)t + G(0) < G{t) < £{0)t^ + F{0)t + G(0); (2.3) 
2m 
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if'y>^, then 

— 12 + F{0)t + G(0) < G{t) < fclllf (0)t2 + F{0)t + G(0). (2.4) 
2m 2 

Proof of Lemma 2.2 Let us note that (1.12) implies £{t) < £{0). 
Identities (2.1) and (2.2) imply for 7 < | 

G"it) = 2Ek{t) + Emit) + 3(7 - l)E,it) = 2Ek{t) + Emit) + 8(7 - l)E,{t) = 

- 2£{t) - E-^{t) - (2 - 3(7 - l))E,it) < 2£(0), (2.5) 

and for 7 > I 

G"it) = 2Ekit) + Emit) + 3(7 - l)E,{t) =. 2Ek{t) + E^it) + 3(7 - l)E,{t) = 

= 3(7 - l)£it) ~ (37 - 4)S„(t) - (37 - h)Ek{t) < 3(7 - 1)£(0). (2.6) 
Inequality (2.5) and (2.6) result in the right-hand sides in (2.3) and (2.4). 

Further, since Ek{t) > |^ (that follows from the Holder inequality), Ei{t) and 
Em{t) are nonnegative, then from (2.1), (2.2) we obtain 

G"it) > — . (2.7) 
m 

Thus, the left-hand sides estimates in (2.5) and (2.6) are obtained. Lemma 2.2 is 
proved. 

The next step is two-sides estimate of Ei{t) + E,n{t). 

Lemma 2.3. If P =/= 0, then for highly decreasing at infinity classical solutions 
to (1.1) - (1-5) for sufficiently large t the following estimates hold: 

G3(.-?)/2(,) ^ W + W ^ G^i.%m) ^ (2.8) 

for 7 < I , and 

^^Sw^<m)+Ernit)<^^, (2.9) 

for 7 > |, with positive constants Ci, G2 (Ci < C2). 

Proof of Lemma ^.5 The lower estimate is due to '8' . For an arbitrary dimension 
of space it follows from the inequality 

^ 2-y ^ "+2)"'"' 

||/||Li(Ri;d2;) < C'7.n|l/llL-'(R";d2;) II/I1l1(R'«'; |a;Pda;)' 

with 

G-y — 

Thus the constant 



27 \ + / 27 



n(7 — 1) / \^(7 ~ 1) 



G, = ^imG-iy-^. 
7 — 1 

The method of the upper estimate of Ei{t)+Em{t) is also similar to [8]. Namely, 
let us consider the function Q(t) = 4G{t)£{t) — F'^it). The Holder inequality gives 

F2 < AG{t)Ek{t), therefore £{t) = Ek{t) + Em{t) + E,{t) > E^{t) + and 

m) + Ern{t)<-^. (2.10) 
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We notice also that Q{t) > provided the pressure does not equal to zero identically. 
Then taking into account (2.1), (2.2) and (1.12) we have 

Q'(t) = AG'{t)£it) - 2G'{t)G"{t) + AG{t)£'{t) = 

= 2 G'{t) ((5 - 37) E,{t) + Em{t)) + AG{t)£'{t) < 

< 2 G'{t) ((5 - 37) E,{t) + E„M- (2.11) 

Since P ^ 0, inequality (2.7) guarantees the positivity of G'{t) for sufSciently 
large t, then for 7 < | 

Q'(i) =2G'(t)((5-37) iE,it)+E^{t)) + iSj-4)E„{t)) < 

< 2 G'{t) (5 - 37) {Ei{t) + E^it)); (2.12) 

and for 7 > | 

Q'{t) ^2G'it)iiE,it)+E^it)) + {A-3j)E^{t)) <2G'{t){Ei{t)+E^{t)). {2.n) 

Thus, for 7 < I inequality (2.12) implies 

Q'{t) 5-3jG'it) 
Q{t) - 2 G{t) ' 

whereas for 7 > | inequality (2.13) results in 

Q{i) - 2 Git) ■ ^^-'^^ 

Therefore (2.10), (2.14) and (2.15) imply the upper estimates in (2.8), (2.9), 

Q(o) 



where C2 = iQ(0)G^(0) in (2.8) and C2 = in (2.9). The proof of 



Lemma 2.3 is over. 



Now we can estimate the rate of decay of total energy £. Prom (1.12), (1.8) we 
have 

£'{t) < -aK{Ei{t))~^(^. (2.16) 
Further, nonnegativity of Em{t), Lemmas 2.2 and 2.3 and inequality (2.16) imply 

£'{t) < -uKG^^^ G^t) <-Lt^, 
with a positive constant i, for 7 < |, and 

£'{t) < -1/ if Ca" G'5l^(t) < -Lt^^(^, 

for 7 > |. In both cases this contradicts to the non- negativity of £{t). However, 
the decay of total energy in the case of highly decreasing solutions is greater that 
follows from inequality (1.13). 

The author thanks Prof. Ronghua Pan for attraction of attention to this prob- 
lem and helpful discussion. 
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